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ABSTRACT 

The  theoretical  model  for  the  dynamic  response  of  multi-directional  functionally  graded  thin 
plates  under  an  in-air  blast  loading  from  a  Friedlander  type  pressure  loading  is  presented.  The 
theory  is  presented  in  the  context  of  the  classical  linear  plate  theory  (CPT)  which  is  based  on  the 
Kirchoff-Love  assumptions.  The  plate  is  assumed  to  be  thin,  in-plane  strains  are  small  compared 
to  unity,  and  the  transverse  and  normal  strains  are  considered  to  be  negligible.  Additionally,  the 
theoretical  model  assumes  that  the  material  properties  of  the  two  constituent  materials  vary  in  all 
three  coordinate  directions.  This  implies  in-plane  as  well  as  through  the  thickness  grading 
according  to  3  independent  power  law  distributions.  Simply  supported  boundary  conditions  are 
assumed  along  all  four  edges.  The  governing  equations  of  motion  are  derived  through  the  use  of 
Hamilton’s  Principle.  The  dynamic  response  is  determined  through  the  use  of  numerical 
integration,  using  the  Gaussian-Quadrature  Method,  the  Galerkin  Method,  and  the  Fourth-Order 
Runge-Kutta  Method  with  zero  initial  conditions.  Results  are  presented  using  the  technique  of 
spatial  tailoring  to  determine  the  optimization  of  the  3D-Grading  from  a  response  standpoint. 
Finally,  validations  are  made  with  simpler  cases  found  within  the  literature. 
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1.  Introduction 

With  the  demand  for  protection  from  the  effects  of  war,  combat  military  vehicles  need  to  be 
fully  armored  to  protect  the  soldier  from  such  things  as  gun  fire  (and  or  penetration),  land  mines, 
IED’s,  Fire,  etc.,  One  of  the  most  daunting  tasks  is  to  engineer  a  military  vehicle  to  withstand 
blast  from  an  IED  or  land  mine.  The  most  important  areas  of  the  vehicle  are  the  hull  which 
resides  underneath  the  vehicle  and  armor  plating  surrounding  various  key  locations  of  the  vehicle 
to  protect  the  occupants  inside.  In  the  case  of  the  hull,  a  structurally  sound  hull  would  protect  the 
occupants  from  large  crushing  floor  accelerations  impacting  the  lower  part  of  the  body.  For  this 
reason,  such  new  materials  as  multidirectional  functionally  graded  materials  should  be  explored 
to  determine  any  possible  advantages  over  other  types  of  materials. 

Multidirectional  materials  are  materials  which  exhibit  properties  which  are  a  function  of  all 
three  coordinate  directions  as  opposed  to  just  the  transverse  direction  found  in  conventional 
functionally  graded  materials.  This  leads  to  the  concept  of  spatial  tailoring  whereby  by  one  can 
manipulate  the  distribution  of  the  constituent  materials  in  all  three  coordinate  directions  to 
achieve  an  optimum  distribution  which  will  provide  the  desired  structural  response.  This  is  the 
focus  of  this  paper. 

Currently,  there  is  a  disparity  of  research  and  information  on  this  topic.  Very  little  exists  in 
the  literature.  Due  to  the  complexity  of  the  theoretical  development,  considering  a  thorough 
treatment  of  the  theoretical  model,  which  most  likely,  would  require  advanced  numerical 
solution  procedures,  a  simplified  model  based  on  the  linearized  theory  is  adopted  to  serve  as  a 
basic  foundation  upon  which  to  build  upon.  Three  independent  power  law  distributions  are 
presented  which  describe  mathematically  the  grading  in  all  three  coordinate  directions.  The 
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transverse  or  through-the-thickness  grading  as  well  as  the  in-plane  grading  is  assumed  to  be 
symmetric  on  both  the  top  and  bottom  faces. 

The  present  linear  theory  includes  damping  effects,  transverse  inertia,  and  a  transient  normal 
loading,  due  to  a  Friedlander-type  pressure-time  impact. 

It  should  be  emphasized  that  this  theoretical  model  is  only  valid  within  the  elastic  region  of  the 
material.  To  remain  in  the  elastic  region,  the  magnitude  of  the  transverse  pressure  acting  on  the 
plate  would  most  likely  need  to  be  of  a  low  to  moderate  intensity.  The  stresses  and  strains  would 
have  to  be  evaluated  to  determine  if  they  reside  beyond  the  elastic  limit. 

2.  Basic  Assumptions  and  Preliminaries 

Shown,  in  Figure  1,  is  a  pictorial  representation  of  a  multidirectional  functionally  graded  plate 
referred  to  an  Orthogonal  Cartesian  Coordinate  System  (x,  y,  z) ,  where  z  is  measured  positive 
in  the  upwards  direction  from  the  mid-surface  of  the  plate.  While,  h  is  the  uniform  thickness  of 
the  plate.  Let  any  two  constituent  materials  comprise  a  functionally  graded  plate.  Then  by 
applying  the  rule  of  mixtures,  a  generic  property  P  ( x ,  y,  z)  can  be  expressed  as 

P(x,  y,  z )  =  P\ U  (x,  y,  z)  +  P2V2  (x,  y,  z)  (1) 

Where  P  represents  the  Young’s  Modulus,  Density,  Poisson’s  Ratio,  Coefficient  of  Thermal 
Expansion,  Etc.,  while,  ,  V2  represent  the  volume  fractions  of  the  two  constituent  materials 
which  must  obey  the  following  relationship 

V{  (x, y,z)  +  V2(x, y, z)  =  1  -  (2) 

With  the  use  of  Eq.  (2),  Eq.  (1)  can  be  expressed  as 

P (x,  .y,  z)  =  [f|  -  P2  \V\  ( x ,  y,  z)  +  P2-  (3) 

Expressing  Eq.  (3)  in  variable  separable  form  gives 
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P(x,  y,  z )  =  [P{  -  P2  ]VCX  (x)Vcy  ( y)Vcz  ( z )  +  P2 


(4) 


Where  the  chosen  functional  relationships  for  the  volume  fractions,  Vcx  (x),  Vcy  (y),  and  Vcz  ( z) 
are  given  in  a  polynomial  and  power  form  as 


r*  ' 

F„U)  =  |  — 

ll< 


1 - 

V 


wN' 


(5a) 


r  f 

I  y  y 

v  (y)  =  I  —  I  -  — 
17  L  L 


(5b) 


Vr7  (z)  = 


/  A 

Z 

1V1 

f  1  +  sgn(  z)7 

r  \ 

-  z 

1V1 

f  1  -  sgn(  z)7 

v  2  j 

|  + 

l  2  J 

(5c) 


The  Signum  function  is  defined  as 


1,  z  >  0 
sgn(  z)  =  \  0,  z  =  0 
-1  z  <  0 


(6) 


N1 ,  N  2 ,  and  M  are  referred  to  as  the  volume  fraction  indexes  providing  a  measure  of  the 
variation  of  the  material  profile  through  the  structure  in  all  three  coordinate  directions.  This 
chosen  grading  of  the  constituent  materials  throughout  the  plate  leads  to  a  symmetric  distribution 
in  all  three  coordinate  directions.  It  should  be  noted,  depending  on  the  3D  grading  desired  that 
other  possible  functional  relationships  for  the  volume  fraction  could  be  utilized. 

The  chosen  constituent  materials  for  this  paper  are  ceramic  and  metal.  This  leads  to  the 
expression  of  the  material  properties  given  by 

[E(x,  y,  z),  p(x,  y,z)]  =  [E cm  ,  p cm  Wi(x,y,z)  +  [E m  ,  pm  ]  (7) 

Where, 
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=  E, 


=  Pc 


(8) 


It  should  be  noted  that  the  variation  of  Poisson’s  ratio,  v(x,  y,  z) ,  is  approximated  as  being 
constant  throughout  the  material  grading  of  the  structure  with  the  assumption  that  the  effect  of 
any  point-wise  variation  on  the  dynamic  structure  response  would  be  minimal  and  or  negligible 
for  a  thin  plate. 

3.  Kinematic  Equations 


3.1  The  Displacement  Field 

Based  on  the  classical  plate  theory,  the  3-D  displacement  relationships  are  expressed  as 

8w  o 

u  -  u  o  -  z - 

dx 

5wq 

U  =  v0-z  — 
dy 

W  =  W  Q 

Where,  (u,  u,  w)  are  the  3-D  displacement  quantities  and  (w0  ,v0 ,  w0 )  are  2-D  displacement 
quantities  of  the  mid-surface  of  the  plate. 

3.2  Linear  Strain-Displacement  Relationships 

Assuming  the  transverse  and  normal  strains  are  negligible,  The  linear  strain  displacement 
relationships  for  plane  stress  are  given  by 
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With  the  use  of  the  displacement  relationships,  Eqs.  (9a-c),  the  strain-displacement  relationships, 
Eqs.  (10a-c,)  can  be  expressed  in  terms  of  2-D  displacement  quantities  as 
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In  the  above  expressions,  (e^.  ,  £^  ,  ) are  known  as  the  membrane  strains;  While, 


(k  ^  ,  k  ,  k  )  are  the  bending  strains. 


4.  Constitutive  Equations 


The  constitutive  equations  for  a  point-wise  isotropic  material  are  given  by 
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Where  the  material  stiffnesses  Qu  (x,  y,  z, ) ,  (i,  j=  1,2,6)  are  given  by 


Sn  =  e22  -  e,2  =  o«  -  a, 6  =02,-0 
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5.  Equations  of  Motion 

Adopting  an  energy  approach,  the  equations  of  motion  are  derived  through  the  use  of  Hamilton’s 
Principle.  It  is  provided  as 


SJ 


-U  +  V)dt  =  0 


(15) 


Where  t0 ,  q  are  two  arbitrary  instants  in  time.  U  denotes  the  strain  energy,  V  denotes  the  work 

done  by  surface  tractions,  edge  loads,  body  forces,  and  damping  forces.  For  this  paper,  there  are 
only  surface  tractions  in  the  form  of  a  transient  transversal  loading  and  damping  forces.  T 
denotes  the  kinetic  energy  of  the  structure,  while  S  is  the  variational  operator.  The  strain  energy, 
U  is  given  by 

h 

5U  =  in  \2h(cTxxSexx  +ayySeyy  +  T  *yr  n  (16) 

2 

Q  denotes  the  mid-surface  area  of  the  plate.  The  work  done  by  the  transverse  pressure  and 
damping  forces  is  given  by 


SV  =  l  [Pt(x,  y,t)  -  ju(x,  y)w 0)]dw 0dQ.  (17) 

Jq. 

In  the  above  expression,  Pt  is  the  distributed  force  at  the  top  surface,  and  u(x,  y)  is  the  damping 
coefficient  per  unit  area  of  the  plate.  It  should  also  be  noted  that  although  n  =  ju{x,  y,  z) ,  the 
damping  coefficient  is  approximated  as  being  constant  through  the  thickness  of  the  thin  plate  but 
can  vary  across  the  2D  plane  of  the  plate.  The  transverse  kinetic  energy  is  given  by 
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ST  =  f  j"  p(x,  y)w qSw QdCldt  (18) 

Jo  Jq 

Where,  p{x,  y ) ,  the  inertia  term  is  given  by 


h 

p(x,y)  =  j2hP(x,y,z)dz  (19) 

2 

where  p(x,  y,  z)  is  the  mass  per  unit  volume. 

With  the  expressions  for  the  strain  energy,  Eq.  (16),  the  work,  Eq.  (17),  and  the  kinetic  energy 
Eq.  (18),  in  hand  and  considering  Eqs.  (11),  (12a, b),  (15),  carrying  out  the  integration  throughout 
the  thickness,  integrating  by  parts  where  ever  feasible,  and  taking  into  consideration  the  arbitrary 
and  independent  character  of  variations  results  in  the  equations  of  motion  and  the  associated 
boundary  conditions  which  are  provided  in  terms  of  stress  resultants  and  stress  couples  as 
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The  associated  boundary  conditions  become: 
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dM  xx  dM  ty 
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The  quantities  with  over-carets  are  prescribed  quantities  along  the  boundary  of  the  plate.  It  is 
desired  to  express  the  governing  equation  of  motion  and  the  associated  boundary  conditions  in 
terms  of  displacements.  To  achieve  this  end,  the  stress  resultants  and  stress  couples  are  defined 
and  expressed  in  terms  of  displacements  as  presented  below. 

The  stress  resultants  and  stress  couple  resultants  are  defined  as 

h 

iN ap’M  *»)=  I\<r4,(hz)dz  ,  (31) 


where  a, (3  =  ( x ,  y) .  Substitution  of  the  constitutive  Eqs.  (13)  and  the  strain-displacement 
relations  Eqs.  (11),  results  in  a  relationship  between  the  stress  resultants  and  stress  couples  in 
terms  of  the  mid-surface  and  bending  strain  components  expressed  as 
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Where,  the  global  stiffness  quantities  are  defined  as 

h 

[Aij(x,y),Bij(x,y),Dij(x,y)]  =  J  2/;  Q tj  (x,  y,  z)(  l,  z,  z2  )dz ,  (i,  j  =  1,  2,  6)  (34) 

2 

It  should  be  noted  that  in  the  case  of  symmetric  grading  in  all  three  coordinate  directions 

[Bjj  ]  =  0  .  As  a  result,  these  quantities  are  excluded  from  Eqs.  (32)  and  (33).  With  the  use  of  Eqs. 

(7)  and  (14),  and  carrying  out  the  indicated  integration  within  Eq.  (34),  the  global  stiffness 
quantities  can  be  expressed  as 


—^(Ei,e2,e3) 

1  -  v 


(35a) 


(35b) 

(35c) 
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(36c) 


Making  use  of  Eqs.  (12a, b),  (32),  (33)  with  Eqs.  (20)-(22)  and  adopting  some  basic  algebraic 
techniques,  yields  the  equations  of  motion  and  boundary  terms,  in  terms  of  displacements  as, 
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Su, 


in  o>  y) 


<,2 

6  Uq 


8x‘ 


+  [An  (x,  y)  +  A66  (x,  y)] 


a2 

g  ^0 

dxdy 


+  A66  (x,y) 


D~uq  dA^(x,y)du  0 


dx‘ 


dx 


dx 


dA66  (x,  y) 


dy 


du  q 
Sy 


dvQ 

dx 


+  SAn  (x,  y)  du0 


dx 


dy 


(37) 


Su  i 


l22 


O,  y) 


d2Vp 
dy 2 


+  [  A12  (x,  y)  +  A66  (x,  y)] 


O  Uq 

dxdy 


+  A66(x,y ) 


d~u0  dA^(x,y)dvQ 


dx ‘ 


dy 


dy 


dA66  (x,y) 


dx 


diiQ  du  o 


dy 


dx 


+  dAn  (x,  y)  du0  _ 


dy 


dx 


(38) 


c  „  ,  >vv4  8Dn(x,y)  5(V2w0)  ^  8Dn{x,y)  0(V2wo)_ 

ow q.  D n  (x ,  y )\  Wq  +  2  h  2  h 


v  2£>u  (x,  y)(  V  2w0)  -  (1  -  v) 


Sx  5x  Sy  dy 

{  d2Pn  (x,  y)  g2w0  02£>n  (x,  y)  02w(- 


Sx‘ 


dy' 


dxdy  dxdy 


d2 D u  (x,  y)  d2w 0  ) 


dy- 


8x‘ 


+  p(x,  y)w0  +  p(x,  y)w0  =  Pt(x,  y,t ) 


with  the  associated  boundary  terms  as, 
Along  the  edges  x  =  (0,  L] ) 


du 


0 


dur 


An  ( x ,  y ) - 1-  A12  (x,  y) - =  0  or  u q  =  z<0 

0x  Sy 


a66  o,  y) 


f 8uq  0do  ^ 
dy  dx 


=  0  or  v0  =  v0 


D  ii  ( x ,  y) 


p2 

O  Wq 

dx2 


d  12  O,  y) 


p2 

C7  Wq 

dy2 


=  0  or 


Sw0  Sw0 


dx  dx 


(39) 


(40) 

(41) 


(42) 


Page  11  of  28 

UNCLASSIFIED:  Distribution  Statement  A.  Approved  for  public  release 


UNCLASSIFIED:  Distribution  Statement  A.  Approved  for  public  release 


8  w  n  8  w  n 

-  Dn(x,  y ) - —  -  [£>12  (x,  y )  +  4 D66  (x,  y)] - -- 

dxi  dxdy 

or  w  =  Wq 

8D ii  (x,  y)  d2w0  /[  8P 66  (x,  y)  d2wQ  dPl2  (x,  y)  d2wQ  _ 
dx  qx 2  5y  dxdy  8x  fiy2 


(43) 


Along  the  edges  y  =  (0,  L2 ) 


A66  A,  j) 


<9m0 


0Uq 

5x 


=  0 


or 


11 0  -  u  o 


0Uq  du  q 

A22(x,y) - +A12(x,v) - =0  or  u0  =  u0 

8y  dx 


D 22 (*. V) 


P2 

o  w  q 

5y2 


£>i2  (x>  y) 


s2 

O  W  0 

5x2 


=  0 


8w0  8wq 
or  - =  - 

8y  8y 


8  Wr,  8  Wft 

-  D22  (x,  y) - —  -  [£>12  (*.  y)  +  4 D66  (x,  y)]  — - - 

8y  dx  8y 

8D22(x,y)  d2w0  8D66(x,y)  d2w0  8DU  (x,  y)  d2  w0 

- - - 4 - =0  or  w0  =  vr  q 

8y  dy 2  dx  dxdy  dy  qx2 


(44) 

(45) 


(46) 


(47) 


.  By  observation,  it  is  seen  that  the  first  two  equations  of  motion,  Eqs.  (20)  and  (21)  governing 
the  in-plane  motion  are  decoupled  from  the  third  equation  of  motion,  Eq.  (22),  which  governs  the 
transverse  or  bending  motion  of  the  plate.  As  a  result  of  the  decoupling,  the  governing  system 
(The  equations  of  motion  and  boundary  terms)  is  reduced  to  one  governing  equation  of  motion 
and  two  boundary  conditions  along  each  edge  of  the  plate  (Note:  a  fourth-order  differential 
equation  requires  two  boundary  conditions  along  each  edge).  This  reduces  the  governing  system 
to 
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4  8Dn(x,y)  0(V2wo)  8Dn{x,y)  0(V2wo) 

owq  :  P\\  (x,  y)V  w q  +  2 - r  2 - h 


V2Dn(x,y)(V2w0)-(l-v) 


0x  dx  dy  dy 

f  g^Pn  P,  y)  g2w0  0~£>n  P’ y)  d‘ 


'0 


dx  ‘ 


dy- 


dxdy  dxdy 


(48) 


d2Pll  (x,  y)  02wo  ] 


dy‘ 


0x" 


+  p(x,  y)w0  +  p{x,  y)w0  =  Pt(x,  y,t ) 


Along  the  edges  x  =  (0,  L{ ) 


~2 

0  W  r 


P2 

0  Wr 


- Dn(x,y ) - r^~£12U,  y) - ^=°  or 

dx  dy 


dw  o  0wo 


dx  dx 


(49) 


P>  y ) 


0  Wq 
0X3 


[Z)12  (x,  y)  +  4 D66  (x,  y)] 


03wo 
0x0y  2 


ggu  (X,  y)  d2w0  4  0£>66  Q’  y)  d2w0  dDn  (x,  y)  02wo  _ 


dx 


8x‘ 


dy  dxdy 


dx 


dy  ‘ 


or  w  =  Wr 


(50) 


Along  the  edges  y  =  (0,  L2) 


p2 

O  XV  ft 

p 22  O.  y) — —  -  D\i  P>  y) 

dy2 


a2 

0  Wr 


dx‘ 


=  0  or 


0Wq  0Wq 


dy  dy 


(51) 


d 22  O’  y) 


n3 

O  XV  q 

dy  3 


[P\2  O’  y)  +  4D66  (x,  y)] 


5  3  Wq 
0x20y 


dP 22  (x,  y)  0  wc 
3y  0y 2 


1  aD66  O’  y)  d2w0  0^12  (x,  y)  02wo  o 


0x  0x0y 


0y 


0x" 


Or  Wq  =  Wq 


(52) 


For  the  problem  at  hand,  simply  supported  boundary  conditions  are  chosen.  This  leads  to  the 
following  choice  of  boundary  constraints  expressed  as: 

Along  the  edges  x  =  (0,  Ll ) 
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wo=0’  M  xx  =  ~Dn(x,y) 


o  w( 


dx' 


-  Dn  (x,  y )  ■ 


O  Wr 


=  0 


(53) 


Along  the  edges  x  =  (0,  L,) 

2  2 

d  w0  d~w () 

w0  =  °’  M  yy  =  ~D22  O’  - --^12^^) - ~  =  0  •  (54) 

dy  ~  dx 

6.  Solution  Methodology 

To  facilitate  a  solution  of  Eq.  (48),  the  Galerkin  Method  is  chosen  which  requires  that  both  the 
essential  (kinematic)  and  natural  boundary  conditions  be  fulfilled.  To  achieve  this  requirement, 
the  transverse  displacement,  w0  (x,  y) ,  is  expressed  in  functional  form  as, 


w0(x,  y,t)  =  W  (r)  sin(  Amx)sin(  /uny)  (55) 

Where  Xm  =  m  x  j L{  ,  jun  =  nx / L2  and  (m  ,n)  are  the  number  of  sine  half-waves  in  the 
corresponding  directions.  A  suitable  representation  for  Pt  ( x ,  y,  t )  is  given  in  terms  of  a  Navier 
representation  as 

Pt{x,  y,t)  =  Pmn  (t ) sin(  Amx)sm(  funy)  (56) 


Integrating  both  sides  of  the  above  expression  over  the  plate  area  gives 


Pmn  (0  = 


4  rL2  rLl 


- f  f  Pt(x,y,t)sm(  A, mx)sm(  junydxdy 

l  L2  Jq  Jq 


(57) 


The  pressure  loading  is  given  by  the  Friedlander  expression  as 


r 


pt(x> y>t)  =  P,(t)  =  (Pso  -  P0)j  i  - 

L 


-c-u/c 


(58) 


Substituting  into  Eq.  (57)  and  integrating  gives 
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Pmn  (t)  =  \6  Pt(t) / mn  k  2  (59) 

In  Eq.  (58),  Pso  is  the  peak  overpressure  above  ambient  pressure,  P0  is  the  ambient  pressure,  ta  is 
the  time  of  arrival,  tp  is  the  positive  phase  duration  of  the  blast  wave,  t  is  the  time,  and  a  is 
known  as  the  decay  parameter  which  is  determined  by  adjustment  to  a  pressure  curve  from  a 
blast  test. 

Applying  the  Galerkin  Method  to  Eqs  (48),  with  the  use  of  Eqs  (55),  (56),  and  (59),  gives 

I  W  ( t)+C  W  ( t)  +  K  w  (t)  =  16 P'^  ,  (60) 

mn  mn  v  7  mn  mn  v  7  mn  mn  v  7  4  7  x  7 

mn  n 


Where  W mn  (. t )  is  the  amplitude  of  deflection,  I mn  is  the  plate  inertia,  C  mn  is  the  overall  damping 
coefficient  of  the  plate,  and  K  mn  is  the  overall  stiffness  of  the  plate  all  of  which  are  given  as 


K, 


Jo  Jo  [fc  +  )  DU  y)sLsny  ~  2A, 


fam  +  Pn  ) 


2  )2  SDn  (*»  y) 


dx 


2 

n  c  c  — 

mx  mx  ny 


(61) 


,2  2  \2dDn(x,y)  2 

1  y 

mx  ny  ny 


-  2Pn  k  +  K  )  —  ",  '  smx  cny  sny  ~  \Kn  +  VP  n  I - + 


fe  +  VP  l ) 


2)5  Dn(x,y)  2  2 

i  "  c 

mx  ny 


dx' 


d~Du{x,y)  l  2  9  \52Z)ii  (x,y)  2  2  1 

+  2 Xmnn  (1  -  v ) -  cmx smx c  sny  -  (//„  +  v4 -  j - smx sny  \ dxdy 

dxdy  -  dy - 


C 


mn 


A  2  U  2  9 

P(x,y)smxs  dxdy 

0  Jo 


(62) 


2  2 

'mn  =  Po(x,y)smxs  dxdy 

Jo  Jo 

Within  the  above  expressions, 


(63) 
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cmx  —  cos(  m  7ix j )  ,  smx  =  sin (m/rx/Lj)  (64a, b) 

cny  =  cos(  nny /l2  )  ,  j ny  =  sin(  nxy/L2)  (65a,b) 

Eq.  (60)  can  be  normalized  by  dividing  through  by  lmn  which  results  in, 

^  mn  (^)  mn  ® inn  wm,At)  + 

® mn^ mn  C'  '  (0  (66) 

where  comn  =  ^Kmn  /lmn  is  the  natural  frequency  of  the  plate,  A  mn  =  Cmn  /2 lmn  comn  is  the 

overall  normalized  damping  coefficient,  and  Pmn  (. t )  =  16  Pt  ( t ) / mn  n 4  /  is  the  normalized 
transverse  pressure. 

7.  Results  and  Discussion 

After  a  thorough  literature  search  was  conducted,  it  was  found  that  there  is  none  to  very  little 
work  commenced  on  multidirectional  functionally  graded  materials  and  or  plates.  This  area  of 
research  is  very  sparse  at  least.  A  couple  of  papers  were  found  that  introduce  the  subject  such  as 
[1]  Birman  and  Byrd  and  [2]  Birman  et  al.  As  a  result,  nothing  can  be  found  in  terms  of  the 
theoretical  developments  to  validate  the  deformation  of  the  plates.  On  the  other  hand,  some 
validations  can  be  made  with  the  natural  frequency  of  the  multidirectional  functionally  graded 
plates  where  the  theory  is  simplified  down  to  simpler  cases.  For  this  reason.  The  governing 
differential  equation,  Eq.  (66),  which  governs  the  structural  dynamic  response  of  multidirectional 
functionally  graded  plates  is  utilized  to  determine  the  natural  frequency  for  the  special  case  of  a 
unidirectional  functionally  graded  plate  at  either  the  fully  metal  and  or  fully  ceramic 
compositions.  For  all  cases  of  comparisons,  Ti-6A1-4V/Aluminum  Oxide  was  used. 

In  Table  1.,  A  comparison  is  made  between  the  present  theory  and  those  reported  by  [3] 
Befarani  et  al.  and  [6]  Bishop,  for  a  simply  supported  square  isotropic  plate,  at  various  modes 
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where  the  frequency  is  non-  dimensionalized.  Very  good  agreement  is  seen.  For  this  isotropic 
case,  (M  =  0, Nl  =  0,  N2  =  0  ).  In  table  2.,  comparisons  are  made  for  the  natural  frequency  in 
Hertz  of  a  unidirectional  functionally  graded  plate.  Comparisons  are  made  for  both  the  fully 
metal  (Ti-6A1-4V)  and  fully  ceramic  (Aluminum  Oxide)  compositions  for  two  modes  of 
vibration  between  the  present  theory  and  those  reported  by  [3]-[6].  These  comparisons  are 
presented  for  the  case  of  fully  metal  (M  =  0 .  N ,  =  0,  A,  =  0)  and  for  the  case  of  fully  ceramic  ( 
M  =  l,Nt  =oo,  N2  =  co  ).  Very  good  agreement  can  be  seen.  The  material  properties  utilized  for 
Ti-6A1-4V/Aluminum  Oxide  are  provided  as, 

Ti-6A1-4V:  E  =  105  .7  x  10  9  Pa ,  v  =  0.298  ,  p  =  4429  kg  /  m  3 

Aluminum  Oxide:  E  =  320  .2  x  10  9  Pa ,  v  =  0.260  ,  p  =  3750  kg  /  m  3 

For  the  present  case,  various  relationships  between  the  geometrical  and  or  material  parameters 
and  their  effect  on  the  structural  response  due  to  blast  have  been  studied.  Fig  2.  Which  compares 
the  central  deflection  as  a  function  of  time  for  the  isotropic  cases  (Ceramic  and  metal)  as  well  as 
for  2  types  of  functionally  graded  materials  (Bi-directional  and  multidirectional).  It  clearly  shows 
that  the  metal  composition  has  higher  amplitudes  and  lower  frequencies  of  oscillations  in  the 
absence  of  damping  than  for  the  cases.  Ceramic  appears  to  be  the  best  performer  from  an 
amplitude  standpoint  but  has  the  highest  frequency  of  oscillation.  The  Multidirectional  and 
Bidirectional  case  fall  in  between  with  multidirectional  being  the  best  performer  between  the 
two. 

Fig  3.  Depicts  how  various  aspect  ratios  effects  the  deformation-time  response  of  a 
multidirectional  functionally  graded  plate  (MDFGP).  The  optimum  response  appears  to  reside 
with  an  aspect  ratio  of  1.  As  the  aspect  ratios  get  smaller  the  deflections  become  larger  with  a 
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lower  frequency.  As  the  aspect  ratios  become  larger  the  deflections  again  appear  to  be  larger  but 
smaller  than  the  case  of  smaller  aspect  ratios.  Also,  the  frequencies  are  lower. 

In  Fig  4.  The  effcet  of  various  amounts  of  damping  on  the  structural  response  provides  results 
as  you  would  expect.  As  the  amount  of  damping  is  increased,  the  oscillations  decay  faster  over 
time.  For  a  fixed  amount  of  damping,  comparisons  on  the  central  deflection  for  the  cases  of  fully 
metal,  fully  ceramic,  bidirectional,  and  multidirectional  are  presented  in  Fig  5.  It  is  apparent 
again  that  like  the  response  in  Fig  2.  That  ceramic  is  the  best  performer.  With  metal  inherently 
being  the  worst.  Again  the  bidirectional  and  multidirectional  cases  are  in  between  with 
multidirectional  superseding  the  bidirectional  behavior  from  a  performance  standpoint  in  regards 
to  the  severity  of  the  amplitude  of  deflection  the  frequency  of  oscillation.  As  previously  seen 
metal  has  higher  amplitudes  with  lower  frequencies  and  vice  versa  for  ceramic  with  all  other 
cases  in  between.  In  the  last  figure,  Fig  6.,  the  effect  of  the  amount  of  decay  of  the  blast  pressure 
has  on  a  MDFGP  is  depicted.  As  the  amount  of  decay  is  increased  it  shows  that  the  response  is 
diminished  somewhat.  To  substantially  diminish  the  amount  of  decay  for  a  MDFGP,  a  much 
larger  decay  parameter  would  have  to  be  used. 

8.  Conclusions 

In  conclusion,  a  linear  theory  of  thin  multidirectional  functionally  graded  plates  has  been 
provided  where  the  grading  occurs  in  all  three  Orthogonal  Cartesian  Coordinate  directions. 
Comparisons  have  been  made  between  plates  of  other  types  of  grading  such  as  isotropic  (metal 
and  ceramic)  and  bidirectional.  It  has  been  found  that  ceramic  is  the  best  performer  under  blast. 
Although  the  frequency  of  oscillation  is  higher  the  deformation  is  greatly  reduced.  Ceramic  due 
to  its  stiffness  id  very  rigid.  On  the  other  hand  depending  on  how  brittle  the  ceramic  is  could  also 
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play  a  factor.  Within  the  functionally  graded  family,  multidirectional  functionally  graded  plates 
seem  to  perform  the  best  as  compared  with  bidirectional  functionally  graded  plates.  Further  study 
should  demand  a  determination  of  the  stress  distribution  within  these  types  of  structures  and 
material  makeup  to  be  the  subject  of  further  research.  It  is  hoped  that  this  current  study  will  lay 
the  ground  work  for  more  complex  theoretical  cases  and  other  types  of  boundary  conditions. 

Most  likely  all  other  cases  will  require  advanced  numerical  techniques  as  part  of  the  solution 
process. 
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Figure  Captions 


Fig  1.  A  pictorial  representation  of  a  simply  supported  multidirectional  functionally  graded  plate 
exposed  to  a  spherical  blast  above  the  plate. 

Fig  2.  Comparison  of  the  amplitude  of  center  deflection  as  a  function  of  time  for  various  types 
of  functionally  graded  plates. 

Fig  3.  The  amplitude  of  deflection  as  a  function  of  time  for  a  multidirectional  functionally 
graded  plate  for  various  aspect  ratios. 
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Fig  4.  The  amplitude  of  deflection  as  a  function  of  time  for  a  multidirectional  functionally 
graded  plate  for  various  amounts  of  damping. 

Fig  5.  Comparison  of  the  amplitude  of  center  deflection  as  a  function  of  time  for  various  types 
of  functionally  graded  plates  for  a  fixed  amount  of  damping. 

Fig  6.  The  amplitude  of  deflection  as  a  function  of  time  for  a  multidirectional  functionally 
graded  plate  for  various  amounts  of  decay  of  the  blast  pressure. 


Tables 


Table  1.  Comparison  of  the  non-dimensional  frequency  m  =  coLr  ^Jl2p^l--v^y/~Eh^  for 

simply  supported  boundary  conditions  of  a  square  plate  (Ll  =  L2  =  0.4  m ,  h  =  0.002  m) 


^ mn 

Present 

Reference  [6] 

Reference^] 

(Dn 

19.71 

19.74 

19.76 

w12 

49.27 

49.35 

49.37 
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0)13  98.55  98.70  98.74 


Table  2.  Comparison  of  the  natural  frequency  co  (Hz)  for  a  simply  supported  square  Ti- 

6A1-4V/  Aluminum  Oxide  functionally  graded  plate  (L1  =  L2  =  0.4  m,  h  =  0.002  m) 


Mode 

Present 

Reference^] 

Reference^] 

Reference^] 

Reference  [4] 

Fully  Metal 

m=l,n=l 

144.89 

143.4 

143.67 

145.04 

144.66 

m=l,n=2 

362.23 

358.42 

360.64 

362.61 

360.53 

Fully  Ceramic 

m=l,n=l 

270.93 

273.906 

268.60 

271.23 

268.92 

m=l,n=2 

677.33 

685.003 

674.38 

678.06 

669.40 

Figures  (In  numerical  Order) 
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Detonation 
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- An=0 
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Time,  t  (sec) 
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Amplitude  of  Center  Deflection,  Wj  j (t)  (cm) 


UNCLASSIFIED:  Distribution  Statement  A.  Approved  for  public  release 


Page  27  of  28 

UNCLASSIFIED:  Distribution  Statement  A.  Approved  for  public  release 


Amplitude  of  Center  Deflection,  Wji(t)  (cm) 


UNCLASSIFIED:  Distribution  Statement  A.  Approved  for  public  release 


